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Abstract 


We study photoproduction during the non-equilibrium stages of the forma¬ 
tion of chiral condensates within the “quench” scenario of the chiral phase 
transition. The dynamics is modeled with a gauged linear sigma model. A 
novel quantum kinetic approach to the description of photoproduction far off 
equilibrium is developed. We find that non-equilibrium spinodal instabilities 
of long wavelength pion fluctuations are responsible for an enhanced photo¬ 
production rate for energies < 80 MeV at order a. These non-equilibrium 
effects lead to a larger contribution than the typical processes in the medium, 
including that of the anomalous neutral pion decay 7r° —> 2y (which is of order 
a 2 ). We follow the evolution of the dynamics throughout the phase transition, 
which in this scenario occurs on a time scale of about 2.5—3 frn/c and integrate 
the photon yield through its evolution. The spectrum of photons produced 
throughout the phase transition is a non- equilibrium one. For thermal initial 
conditions at the time of the quench it interpolates between a thermal distri¬ 
bution about 6% above the initial temperature (at the time of the quench) for 
low energy < 80 MeV photons, and a high energy tail in thermal equilibrium 
at the initial temperature, with a smooth crossover at 100 MeV. The rate 
displays a peak at ~ 35 MeV which receives a larger enhancement the closer 
the initial temperature at the time of the quench is to the critical tempera¬ 
ture. It is found that the enhancement of photoproduction at low energies is 
not an artifact caused by the initial distribution of the photons, but is due 
to the pionic instabilities. We suggest that these strong out of equilibrium 
effects may provide experimental signatures for the formation and relaxation 
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of DCC’s in heavy ion collisions. 
ll.30.Qc, 11.30.Rd, ll.10.Ef 


Typeset using REVTjrX 


2 



I. INTRODUCTION 


Considerable interest has been sparked recently by the possibility that disoriented chiral 
condensates (DCC’s) conld form during the early stages of evolution in high-energy, high 
luminosity hadron or heavy-ion collisions since this would be a signal of the chiral phase 
transition HI- The main idea is that near the central rapidity region in these collisions, 
typically a large amount of energy ~ a few Gev is deposited in a volume of order a fm 3 , 
corresponding to temperatures above 200 MeV at which the chiral symmetry is restored 
m- It is conceivable that as these regions cool down through hydrodynamic expansion 
they might give rise to domains where the order parameter of QCD points in a different 
direction from the a one, resulting in a disoriented chiral condensate (DCC’s) [l|-f|]. It has 
been speculated that these regions will then act as strong pion lasers, relaxing to the true 
physical vacuum on time scales of a few fm/c by emitting coherent bursts of pions with 
very definite isospin correlations This phenomenon would be a striking experimental 

signature of the chiral phase transition and could provide an explanation for the Centauro 
and anti-Centauro (JACEE) cosmic ray events |5|.|6j. [T0| . |Tl| . 

Since trying to understand such phenomena by investigating QCD is currently impossible, 
efforts have mainly been focussed on the dynamics of low energy effective theories of hadrons, 
some popular examples being the linear and non-linear a models. The rationale for using 
a low energy effective theory approach is that, as in the study of critical phenomena, one 
expects that if a large amplitude coherent configuration of soft modes is produced during 
the transition, such an occurence should be universal and only very weakly dependent on 
the detailed features of the theory. 

Rajagopal and Wilczek JT2[] have argued that the 0(4) linear a model encodes the relevant 
low energy chiral phenomenology of QCD. In particular they argue that the model lies in 
the same static universality class as QCD with 2 light flavours, and that JT3| large coherent 
regions of correlated pions can form only if the phase transition is strongly out of equilibrium. 
They have proposed the “quench” scenario in which the initial “hot region” supercools 
through the phase transition via hydrodynamical expansion on time scales much shorter 
than the relaxation times for long-wavelength fluctuations. As a result, the system enters 
into the spinodal region which induces the growth of the unstable long-wavelength modes 
just as in the process of spinodal decomposition in phase separation 0- This version 
of a disoriented chiral condensate is somewhat different from the “baked alaska” scenario 
envisaged by Bjorken,Kowalski and Taylor (3|-|5| although the relaxation process of both 
configurations is similar. 

Results obtained from studies of both the classical [ 13, 1^-17] and quantum linear sigma 


model [ |TS| , |HJ| , including hydrodynamical expansion ||TB|] seem to suggest that there is a 
range of initial conditions for which there is amplification of pion fluctuations that could 
lead to the formation of large coherent domains during the chiral phase transition. However, 
pions are strongly interacting hadrons, with a typical mean free path in medium of about 
2-3 fm [Q under conditions to be achieved at RHIC or LHC. This results in strong final 
state interactions, and it could be that a DCC signal may be indistinguishable from the 
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background. 

There are currently two experiments searching for a DCC signal, T-864 test /experiment 
at the Fermilab Tevatron collider |21[ and WA98 at CERN-SPS [22j and a proposal for a full 
acceptance detector at LHC (FELIX) [^2[. It is thus imperative to identify an alternative 
signature of the chiral phase transition capable of providing information on the early time 
evolution of the collision and the formation of DCC’s or any other out of equilibrium effect 
associated with the transition. Photons and dileptons are the most promising probes of the 
early stages of the transition. Since they interact mainly electromagnetically, with mean free 
paths typically much larger than the size of the fireball they carry undistorted information 
of the initial non-equilibrium stages of the evolution of the plasma. 

In this context, it is important to realize that there are many sources of photons in rel¬ 
ativistic heavy-ion collisions: (i) direct QCD (hard) photons originating from hard partonic 
processes (typically Drcll-Yan, and qq —> gj), (ii) thermal and non-thermal photons orig¬ 
inating from the quark-gluon plasma, (iii) photons originating from the hadronic gas [^3 


It is the last process that is of interest to us, and is usually the dominant process in the 
production of long wavelength photons with energies less than a few hundreds of MeV’s. 
Some of the main sources of photons in this energy range are the decay of the neutral pion 
n 0 —> 77 , g and rf and resonance decay into photons. Of these, the leading decay channel 
for low energy photons is the “anomalous” decay of the neutral pion. 

In this article we concentrate solely on photo-production, with the particular goal of 
obtaining an ab initio , real-time expression for the photoproduction rate in strongly out of 
equilibrium situations. Our main premise is that if the chiral phase transition is “quenched” 
as envisaged in the scenario of Rajagopal and Wilczek, with the long-wavelength pion modes 
becoming unstable, they will grow almost exponentially in the initial stages of the transition, 
so that these strongly out of equilibrium processes must result in enhanced photoproduction 
at low energies. 

We note that although we focus on the “quench” scenario within a specific phenomeno¬ 
logical model, the new kinetic approach developed in this article is more general and not 
tied to this scenario. In particular, if the transition is strongly first order and is driven by 
nucleation and growth of hadronic bubbles, as suggested by Kapusta and Vischer [^4]], this 
non-equilibrium process should also lead to enhanced photoproduction. 

Dilepton and photon production from a classical coherent pion source have been studied 
by Huang and Wang |25| . These authors studied photoproduction due to a classical Blaizot- 
Krzywicki configuration and found a strong dependence on the initial conditions. However, 
to our knowledge there has not yet been an attempt to compute the photoproduction rate 
including quantum and thermal effects and strong non-equilibrium effects as should result 
from the fast relaxation of DCC’s. It is the point of this article that such a study requires a 
novel formulation of the kinetics of photoproduction (and eventually of dilepton production). 

Typical calculations of photoproduction rates assume a system which is in equilibrium or 
quasi-equilibrium (steady state), and show that the rates are related to the Fourier transform 
of the hadronic current-current correlator pd| [IT . 


The main idea behind our approach is the calculation of expectation values of time 
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evolving observables such as the photon number from first principles, using real-time non¬ 
equilibrium quantum field theory |34|| treating the pion instabilities non-perturbatively 
& These non-perturbative methods were introduced within the context of DCC formation 
in references ||T8l.[TTi|| . Such framework allows us to go beyond the standard results for the 
emission rates to study strongly out of equilibrium situations. In particular we find that 
the emission rates are given by the hadronic (pion) current-current correlation functions, 
but these are obtained from the full non-equilibrium dynamics of the pion fields. Thus the 
emission rates and consequently the photon spectrum obtained using this approach reflect 
all the long range correlations that were generated during the phase transition. Moreover, 
the rate equation resembles a generalized kinetic equation that encodes all the non¬ 
equilibrium time scales of DCC formation, growth and relaxation. In particular, we notice 
the absence of energy conservation at short time scales ~ 1 fm/c, and hence photoproduction 
proceeds via processes that would be kinematically forbidden in an equilibrium situation. 

The focus of our article is to study the photoproduction rates from the coherent DCC 
regions which are formed during the non-equilibrium stages during a ‘quenched’ chiral phase 
transition including the quantum mechanical and thermal effects. 

In section II we give a summary of the non-equilibrium dynamics of the linear sigma 
model. We do not incorporate expansion at this stage, and treat the dynamics of the 
transition via a “quench”. In section III we “gauge” the sigma model and discuss its range 
of validity. In section IV, we present the details of our novel approach to the kinetics of 
photoproduction as an initial condition problem. We obtain the remarkable result that 
in a situation far from equilibrium, as is the case under consideration, the rate is O(a) 
(with the full quantum mechanical current-current correlation function). In section V we 
discuss the numerical results of the calculations of photoproduction, and show that the 
distribution function of produced photons is out of equilibrium interpolating between two 
local equilibrium situations at different temperatures, with an enhanced photoproduction 
rate for energies below 80 MeV. We also compare the results to the case in which there 
are no photons initially and conclude that the non-equilibrium enhancement has its origin 
solely on the pionic instabilities and is insensitive to the initial photon distribution. The 
resulting photon spectrum is out of equilibrium and non-thermal. We compare our results to 
those obtained recently on photoproduction from the hadronic gas |3^j3l|| . Such a comparison 
unequivocally reveals the non-equilibrium features associated with the dynamics of the phase 
transition. In Section VI we discuss the contribution from 7r° —> 2y out of equilibrium and 
argue that it is subleading compared to the lowest order contribution from the charged 
pions. In section VII we present our conclusions, summarize the main results and pose new 
questions. Two appendices are devoted to technical details and a pedagogical example of 
the novel kinetic approach. 


II. THE 0(4)LINEAR SIGMA MODEL OUT OF EQUILIBRIUM 

The linear sigma model has been one of the most popular of the phenomenological models 
to describe the chiral phase transition because it incorporates the correct symmetries and 
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current algebra relations. Obviously such a simple description cannot capture all of the 
details of the underlying theory, QCD, but if the chiral phase transition entails the formation 
of long-wavelength, coherent pion clouds, one expects this to be a robust feature, fairly 
independent of the detailed confining dynamics. 

This model has been studied intensely from the classical, ||T^ . |To| . [To| . |T7[j semiclassical [fT5 |, 
and quantum perspectives both with [|18|] and without [[H| expansion. In this article we do 
not incorporate either longitudinal or spherical expansion |l6|jr8|| , but prefer to illustrate the 
main concepts in the simplest of settings, that of a “quenched” transition. If (as we will see 
later) interesting new physics emerges, there would be motivation to include expansion and 
other features such as inhomogeneities in the analysis of photoproduction from the dynamics 
of the Disoriented Chiral Condensates. Here we summarize the main features of the dynamics 
in the “quench” scenario of the ungauged linear sigma model, and postpone until the next 
section the subtle details of its gauging. 

The Lagrangian density for the linear sigma model with an explicit symmetry breaking 
term is given by 


C - ■ 9*$ - -m 2 (t)§ - + A(l> ■ $) 2 - ha 


( 2 . 1 ) 


where $ is an 0(N + 1) vector, <f> = (cr, if) and tt represents the N pions, tt = 
(7T°, 7T 1 ,7T 2 ,7T 3 , ....tt^ -1 ), and m 2 {t ) introduces phenomenologically “by hand” a quench situ¬ 
ation. Comparison between the “quench” approximation |19| and more realistic treatments 


including longitudinal and spherical hydrodynamical expansions reveal |T8|Jl6|] that there is 
no great discrepancy between the different results insofar as the time scales and correlations. 
We then take 


m 2 (t ) = 


Ml 


rp2 

± c 


0(-f) 


;Ti>T c 


( 2 . 2 ) 


corresponding to a “quench” from the high temperature phase above the critical temperature, 
to zero temperature. 

Although we will be dealing with the 0(4) model with 3 pions, we will use the large N 
limit in order to provide a consistent, non-perturbative framework to study the dynamics. It 
can be implemented by use of an auxiliary field consistently to all orders in 1/iV |35|]. This 
approximation is renormalizable, energy conserving and maintains all of the Ward identities 
associated with chiral symmetry and current algebra. 

As a result, the pions in the asymptotic equilibrium state are exactly massless in the 
absence of the explicit symmetry breaking term. A consequence of Goldstone’s theorem is 
that the spinodal line coincides with the coexistence line, so that there are long-wavelength 
instabilities for all values of the order parameter between zero and the equilibrium minimum 

U- 

The non-equilibrium equations of motion are obtained via the tadpole method 
shift a by its expectation value in the non-equilibrium state 


37 : first 


a(x, t ) = (j>{t) + y(£, t) ; <j>{t) = {a(x, t)) 


(2.3) 
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and then implement the tadpole condition 


(x(x,t)) = 0 ; (9(x,t)) = 0 


(2.4) 


to all orders in the corresponding expansion. The expectation values of the Heisenberg 
operators %(x, t) and tt(x, t ) above are to be computed in the initial density matrix, or alter¬ 
natively, we can compute these expectation values using the Schrodinger picture operators 
and the time-evolved density matrix, as described below. 

To leading order in 1/N the auxiliary field method 
Hartree like factorization 0: 


351 is equivalent to the following 


X 4 -> 6(X 2 )X 2 + constant ; y 3 -> 3(y 2 )y 


(2.5) 


(if ■ if) 2 —> (2 + —){'k 2 ) ; k 2 + constant 


( 2 . 6 ) 


*Y - 7f 2 C\' 2 ) + YY ; A: - (7f 2 )x. (2.7) 

Non-equilibrium quantum field theory requires a path integral representation along a 
complex countour in time |3^j3^| , with the Lagrangian density along this contour given by 

C neq = £[$+] -£[$"] (2.8) 

with the fields < h ± defined along the forward (+) and backward (—) time branches. For further 
details see references [[35|,[!7]]. In the leading order in 1 /N, the effective non-equilibrium 
Lagrangian density is given by 

£[0 + Y + ,7f + ] ~JZ[(f) + x~, 7T _ ] = {^(dx + ) 2 + ^(dn + ) 2 

-V\t)x + - \M x +2 (t)x +2 - \m Y(f)vf +2 } - {+ — -} ■ (2.9) 

where 

V'(t) = 4>{t ) + 0(t)[m 2 (t) + 4A0 2 (t) + 4A(7r 2 )(f)] — h (2.10) 


M%{t) = m 2 (t) + 4A 4> 2 {t) + 4A<7f 2 ) (t) 


( 2 . 11 ) 


M\(t) = m 2 (t) + 12A0 2 (t) + 4A(7f 2 )(f). (2-12) 

Comparing eqns.( pTC| , |2.11|) when 0 = 0 (the equilibrium case) one finds the PCAC 
relation 


/,«; = h 


(2.13) 
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with f n the equilibrium expectation value of the sigma held and M n the equilibrium pion 
mass. The validity of this relation in this approximation is a consequence of the large N 
expansion. A naive Hartree factorization would violate this fundamental Ward identity. 
To leading order in the 1/N expansion, the theory is quadratic at the expense of a self- 
consistent condition and the fluctuations of the sigma held do not influence the dynamics 
of the expectation value or the pion helds and thus decouple from the dynamics. The pion 
holds obey a linear Heisenberg evolution equation (in terms of the self-consistent held), and 
the Heisenberg held operators can be expanded as 

f (f, t) = * + a{Ul(t)e-* *] (2.14) 

vH k ^2 W k ,i 

where a*,, a\. are the destruction and creation operators of Fock states associated with the pion 
held, is the quantization volume, and the mode functions U k (t) and the order parameter 
obey: 

0(f) + [m 2 (f) + 4A0 2 (t) + 4A(7f 2 )(t)]0(t) — h = 0 ; 0(0) = 0; 0(0) = 0 (2.15) 


^ + k 2 + m 2 {t) + 4A0 2 (t) + 4A(7f 2 )(t) 

Uk{ 0) = 1; t4(0) = —iWk,i ; W k) i = 


74(f) = 0 
k 2 + m 2 {t < 0) 


(2.16) 

(2.17) 


(* 2 m=N j 


d 3 k |74(f)| 2 -l 
(2vr) 3 2 W kti 


coth 


W k , 


2 T, 


(2.18) 


Our initial conditions on the pion modes are that at the time of the “quench”, the 
pion modes were in local thermodynamic equilibrium at the initial temperature T, with the 
effective mass m(t < 0) given by ( |2.2| ). The initial condition on 0 is chosen so that the 
expectation value of the sigma held was at the “top” of the potential hill, and the self- 
consistent fluctuation has been subtracted at the initial time , thus renormalizing the mass 
(for more details see reference ||19|| ). 

The linear sigma model is a phenomenological model, and its parameters are hxed by the 
low energy pion physics to be: 


M a ~ 600 MeV ; ~ 95 MeV ; A ~ 4.5 

h « (120 MeV ) 3 ; T c « 200 MeV . (2.19) 

The critical temperature T c is a consequence of the parameters of the model and is 
somewhat larger than the lattice estimates T c ps 150 MeV. The sigma model must be treated 
as a cutoff theory with an ultraviolet cutoff of the order of A ps 1 — 2 GeV. There are two 
reasons for this cutoff: i) for such a large value of the coupling constant the Landau pole is 
at the order of this scale if the parameters (the coupling) are determined at the scale ~ M a , 
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ii) more fundamentally, the simple linear sigma model does not incorporate higher mass 
hadrons, such as vector meson resonances (with masses > 770Mev) and nucleons. The time 
evolution will only be sensitive to the scale of this cutoff for time scales much smaller than 
about 0.1 fm/c, but we are interested in longer time scales of a few frn/c and the numerical 
results on these time scales proved to be insensitive to the values of the cutoff in the range 
1-2 GeV. 

The 

self-consistent numerical solution to the coupled non-linear equations (|2.15| , [2.1^ . p.l7| , ^T8|) 
with the above values of the parameters leads to a complete description of the dynamics. 
The non-equilibrium Green’s functions for the pion fields, which will be necessary in the cal¬ 
culation of photoproduction, are completely determined by the mode functions U k (t) which 
are solutions to the set of self-consistent equations ( |2.16| , [2.1i*| ). In particular, the spatial 
Fourier transform of the necessary Green’s functions are given by : 


(vr+ (k, t)i r+ (-k, t')) = -iS a>b [G£ (t; t')0(t - t') + G£(t; - t)] 


( 2 . 20 ) 


(t T a (k,t)7r b ( -k,t')) = -iS a , b [G£(t;t')0(t' -t) + G£(t; t')Q(t - if)} 


( 2 . 21 ) 


(tt+ ( \k,t)n b {+ \-k,t')) = -iSa }b Gk {>) 


<(>)/+. 


( 2 . 22 ) 


As a consequence of assuming that the system was initially in local thermodynamic equilib¬ 
rium at the temperature T t — 1//3, the Green’s functions G < ,G > obey the KMS condition 


Gk(to] t ) — Gfc (t 0 — ifi', t ). 


(2.23) 


The fundamental Wightman functions are constructed from the pion mode functions, 

solutions of (|2.16|) as: 


G kM = ^-[(1 + n b )U k m*k(t , )+n b U* k (t)U k (t')] 


(2.24) 


G<(M') = --[(1 + n b )U* k (t)U k (t’) +n b U k (t)U* k (t’)] 


2W k 


n b = 


,PWk,i 


(2.25) 

(2.26) 


The numerical results of the integration of the self-consistent evolution equations given 


above have already been discussed in the literature [|19|]. The main feature of the dynamics 
is that for early times t > 0 a band of long wavevectors “see” an inverted harmonic oscillator 
and grow almost exponentially until the backreaction of the fluctuations and the evolution 
of the order parameter shut these instabilities off. These are the “spinodal” instabilities 
that drive the phase transition out of equilibrium. The numerical analysis reveals that 
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at early times long- wavelength pion inodes with wavevectors k < 200 MeV grow almost 
exponentially, whereas the zero mode of the sigma field, i.e. the order parameter, rolls 
towards its equilibrium value close to f n on a time scale ~ 3 fm/c. Figure 1 shows the 
evolution of /(r) = (p(t)/MF vs r = cMpt (Mp = 200 MeV = 1/fm). As shown in figure 1, 
at time scales of the order of 3 fm/c the order parameter reaches f n and begins to oscillate 
around this value, indicating that the phase transition is almost complete. Parametrizing 
[4(f) = \Uk{t)\e~ llfk ^\ we show in figure 2.a ln(|[4(t)|) vs r (time in units of fm/c) for 
several values of k. Figure 2.b shows the values of ^(f) vs r for the same values of k as 
in figure 2.a. Clearly there is a band of long-wavelength pion modes k < 150 — 200 MeV 
that grow almost exponentially during the time for which the order parameter rolls down 
the potential hill towards f n and their phases vary smoothly. The reason that we focus on 
the phases will become clear when we derive the kinetic equations for photoproduction. In 
short, since the phases vary rather smoothly and monotonically during this interval, there is 
no large phase cancellation due to “dephasing” of these modes so that they will contribute 
coherently to photoproduction. These long wavelengths unstable modes do not “decohere” 
during the relevant time scales and as we will see when we study the kinetics, they give the 
primary contribution to the non-equilibrium photoproduction rate. 


III. GAUGING THE LINEAR SIGMA MODEL 


Both charged and neutral pions couple to electromagnetism. Whereas the charged pions 
couple through the charge form factor, the neutral pions couple through the chiral anomaly. 
In the vector meson dominance approximation, the charge form factor of the (charged) pion 
is determined by the exchange of a p vector meson and given by |38| 


G n (q 2 


q 

1 H- y A 

rrip 


(3.1) 


with m p ~ 770 MeV is the mass of the p vector meson. 

Since we will be focusing on on-shell photoproduction of low momentum photons we can 
approximate the charge form factor by 1 to the order that we are calculating (no off-shell 
photons in intermediate states). 

The coupling of the neutral pion is determined by the anomalous quark triangle diagram, 
which for three colors and two quark flavors (up and down) gives the effective “anomalous” 
vertex 


Ca = ^re^F^Fap 7T° 

The physical charged pion fields are given by the combinations 

7r 1 + in 2 7T 1 — in 2 


7T+ = 


y/2 


7T_ = 


V2 


(3.2) 


(3.3) 
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whereas the neutral pion is identified with the third isospin component 7T° and the effective 
lagrangian with electromagnetic coupling in this approximation is given by 


C = + j(*! + *l)A„A» + - jf^F'" + £.4 - V(<r, if) (3.4) 

J/4 = e(7Ti9 M 7r 2 - TTs^TTi) (3.5) 

with V(a,7 r) is given by the potential term in ( |2.1|) in terms of a, 7 ^ 2 , vt°. The current j /t is 
identified with the third component of the isospin current. 

After performing the shift a —> 0(f) + y(T,f), the total effective non-equilibrium La¬ 
grangian is given by 

Cneq = £[0 + x + , ff + , dl+] - £[0 + X", 7f“, = ^(V) 2 + ^(<9tT + ) 2 (3.6) 

- x + 0 - \f+^F+ + e -(7Tt 2 + 7T 2 +2 )A+A + ^ + j+A + » + £ A (vr 3 + , A+) 

- h / ((T,f + )x + + ^y"(<T,7f + )x +2 + ^f /[3l (^ 7?+ )x +3 + ^V [4] (a,n + )x +4 


The fields 7T + should not be confused with the charged pion 7T + ; the superscript in the 
above expression refers to the time branches. However before performing calculations with 
the above non-equilibrium Lagrangian we must address the issue of gauge invariance and 
determine the physical observables. 

Since we want to avoid potential ambiguities with gauge fixing terms and gauge choices, 
we will describe our quantization procedure solely in terms of physical degrees of freedom. 
This is best achieved by passing to the Hamiltonian formulation, using it to recognize the 
physical degrees of freedom and then casting the non-equilibrium calculation solely in terms 
of these, thus avoiding any potential ambiguity with gauge artifacts. 


A. The Electromagnetic Sector: 


Since we are dealing with an Abelian gauge theory, there are only two first class con¬ 
straints for the gauge sector: vanishing canonical momentum for Aq and Gauss’s law (the 
generator of time independent gauge transformations). We can work in terms of gauge in¬ 
variant observables |ID| by projecting on the states in the Hilbert space that are annihilated 
by these first class constraints, or alternatively fix Coulomb gauge which is a physical gauge 
displaying the two transverse physical photon polarizations and the instantaneous Coulomb 
interaction. 

Our first goal is to express the number operator for asymptotic, physical photons with 
two massless transverse degrees of freedom, in a convenient form that can be inserted into the 
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non-equilibrium path integral. In a plasma the “longitudinal” component (the instantaneous 
Coulomb interaction) is screened by a Debye screening length or electric mass, but the 
transverse components are not screened as there is no magnetic mass in the Abelian theory. 

We concentrate only on the physical components i.e. the transverse components of the 
photon held A T (x,t). If II T represents the momentum conjugate to A T (x,t), we can write 
the Hamiltonian for the free electromagnetic theory as follows: 


H = 



-n T • n r + -(v x A T y 


(3.7) 


We then define the Fourier components of the fields as 


P T (k) = 


d 3 x 


n T (F)e 


—ik’X 


(3.8) 


(k) = 


d 3 x 


A T {x)e 


—ik’X 


(3.9) 


where fl is the quantization volume. In terms of these variables the Hamiltonian becomes 

1 


H = j £ [Fr(fc) ■ Pr(-k) + k 2 $ T (k) ■ $ T (-k) 


(3.10) 


Thus the number operator for photons with momentum k per polarization is given by the 
average over the two polarizations 




l 

4 k 


P T (k) ■ P T {-k) + k?$ T (k) ■ $ T (-k) 


1 

2 ' 


(3.11) 


This expression, which is more amenable for use in the non-equilibrium formulation is equiv¬ 
alent to the familiar one given by 




1 

2 


E^(A) 




(3.12) 


which gives the average number of transverse photons per polarization. Here A indicates the 
transverse polarization states, while bt and are the creation and destruction operators 

for photons of momentum k and polarization A. The equivalence is made manifest through 
the mode expansions: 


A T (x,t) 


1 e\(k) 

Vnj- x V2k 



e **+ 6 I-(A) 



(3.13) 


Pt(x, t ) 






k, A 




(3.14) 
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Where the e’s are the transverse polarization vectors. Combining equations (0), 0, 
( |3.13| ), (|3.14|) , we obtain the Fourier components of the transverse phase space variables to 
be 


a L 


4a>W) + 


(3.15) 


P T (k) = b f.(xMk) - 


(3.16) 


These identihcations lead at once to the equivalence between the number operators given by 
( ])• 1 H - fUl up to a normal ordering constant. 

In order to obtain an equation of motion for this number operator we now must address 
the interactions. 


B. The Charged Sector of the Sigma-Model 


In this section we will focus on the interaction between the charged pions and the electro¬ 
magnetic held, under the assumption of a local interaction (charge form factor approximated 
by one). We will postpone the discussion of the anomalous coupling of the neutral pion to a 
later section, where we will argue that such a contribution is subleading out of equilibrium. 
The Hamiltonian (in the gauge invariant sector) for the charged sector is given by 


H — J d 3 x ^n T • fir + ^(V x At) 2 + (Vtt + ) • (Vtt") + e 2 A T ■ A T tt+tC 

+ieAx ■ |7r + V7r - — 7r~V7r + l + Coulomb term 


(3.17) 


The Coulomb interaction will appear in the photoproduction rate at 0(a 2 ) but not to 
order (a); we will neglect this term in our lowest order calculation. A Hartree factorization 
of the “seagull term”, consistent with the large N approximation, leads to a time dependent 
“mass” term for the transverse components 


e 2 A T ■ A t n + n = ^(e 2 A^(r 1 + irf)) ^ l(e 2 ^.(uf + 7r|>) 

= (e 2 4«» = At )4, 


(3.18) 

(3.19) 


The effect of this time dependent “mass term” is to “squeeze” the quantum state for 
the transverse photon states, contributing to photoproduction. However, as it will be clear 
later, such a contribution is of 0(a 2 ) and thus higher order in the electromagnetic coupling 
as compared to the lowest order process that will be seen to dominate the non-equilibrium 
photoproduction rate. 

The electromagnetic current can be rewritten in terms of the pion fields and 7 t 2 using 

0 : 
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j(x,t) = e[f 2 Wi - 7riV7r 2 ]. 


(3.20) 


In terms of the spatial Fourier transform of the pion and gauge fields, we obtain the following 
interaction vertex in momentum space 


J cl 3 xj.A T = \~p) • ; J\-p) 

p 


2 ie 

IT 


^7Ti(M)7r 2 (-£ 

p,k 


p,t){k-$ T (p,t)), (3.21) 


where we have used the fact that p- ^(p, f) — 0. The interaction Hamiltonian, including 
the time dependent mass term for the photon arising from the large N factorization of the 
seagull term, but neglecting the Coulomb interaction is now written as 


H = J2 


^P T {q) ■ Pri-q) + ^{t)$ T (q) ■ ^r(-q) + J(~q) ■ $r(g) 




(3.22) 

(3.23) 


leading to the Heisenberg equations of motion: 


®T(q,t) = P T (q,t) (3.24) 

Pr(q, t ) = -Uf{t)$ T (q, t ) - Jr(q, t). (3.25) 


Here, Jt denotes the transverse component of the current, obtained by projecting the full 
current onto the transverse polarization states: 


J Ti (k) = V il (k)J l (k) ; V«(k) = 6 u 


kih 
k 2 ' 


(3.26) 


Using the definition of the number operator (|3.11|) and the Heisenberg equations of motion, 
we obtain 


Nq(t) — — — p 2 {t){$T{q) ■ $r( — q) + &t(<T) • $t(— q)) + ^r(q) • Jr{—q) + Jr(q) • < £t(— 

4g L 

(3.27) 


The expectation value of this Heisenberg operator equation in the initial density matrix can 
be written in a compact notation as 


(NM = - 


]__d_ 

2 qdP L 
p 2 (t) 9 


j t=t n 


2 q dt" L 


(J^ (g,t) • $ T ( )) 

($r(<fU) ' (-&*") + ($£(&*") ' (-&*)) 


t=t" 


(3.28) 


where we have used fields and currents defined on different branches (forward and backward) 
to avoid the potential ambiguities associated with Schwinger terms in the time derivatives 
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of correlation functions. The expectation value of the emission rate can now be obtained 
by performing a perturbative expansion in a of the expectation values above in the non- 
equilibrium generating functional 


(Nf) = J V[$$}V[$T]V[n+]V[n-]e i{s °- s ° ] e ils t- s i ] N^t) 


with 


st = - my : 




$t(“P) ' $ t(p) + Jt{~p) ■ $t(p) 


(3.29) 


(3.30) 


on the respective ± time branches. The Coulomb interaction can be incorporated in the non¬ 
equilibrium path integral by introducing an auxiliary field variable to replace the Coulomb 
interaction by a linear coupling of the charge density to the auxiliary field (without kinetic 
term), and carrying out the path integral over the canonical momenta of the pion fields. 
However we omit this term since its contribution will be of higher order in a. 

Perturbation theory in a is now carried out in terms of the following Green’s functions 
for the photon field 


<$“*(£, t)$ +J '(-fc, t')) 
($ +i (k,t)$- j (-k, t')) 


^(k)g>M 

v ij (k)^r [ e -^-*')(i + N k ( 0)) + e ik ^N k ( 0) 

z/c 

-iV^{k)G<{t,t') 

V ij (k)— + N k m + e- ik( t- e) N k ( 0) 

2k l 


<$ +i (M)* +J '(-M')> = -iP ij (k) - t') + - t)] 


(3.31) 

(3.32) 

(3.33) 


($- i (M)<i>^(-M')> 
N k ( 0) 


~iV ij {k) [G k (t,t')<d(t - t') + G k {t,t')Q{t' - t)] 

1 

eP k - 1 ’ 


(3.34) 

(3.35) 


where we assume that the initial state is that of local thermodynamic equilibrium and N k ( 0) 
is the initial photon distribution at time t = 0. We will consider two cases: i) the ini¬ 
tial distribution is that of equilibrium at an initial temperature T % = 1//3 and given by 
the equilibrium Bose-Einstein distribution functions, ii) no photons in the initial state, i.e. 
N q ( 0) = 0. We choose this latter initial state for comparison purposes. The case of any 
arbitrary initial distribution can be treated similarly. 

To zeroth order in e, the rate vanishes, as can be seen by taking the time derivatives of 
the free-held propagators ( |3.31| - j3T34| ) . The lowest order contribution to the last two terms in 
( |3.2fj| ) comes from the seagull term which is O(a). Thus, the term proportional to /i 2 (t) will 
give a contribution of order a 2 , since /i 2 (t) oc a, while the current vertex will yield the O(a) 
leading contribution to the photoproduction rate as shown below. This means that we need 
only calculate the current vertex contribution. 
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The first term in (|3.28| ) receives contributions from the J ■ At term in the action. In 
particular, these contributions are proportional to the hadronic (pion) current-current cor¬ 
relators which are 0(e 2 ) and are therefore the leading contributions to the photoproduction 
rate. It should be noted, however, that in an equilibrium situation the O(a) (lowest or¬ 
der) contribution from the current-current correlator vanishes due to kinematics as it will 
be shown explicitly below, and the lowest order contribution to the rate is from the seagull 
terms along with the Compton scattering and pair annihilation diagrams. 


IV. PHOTOPRODUCTION 

As mentioned above, photons and dileptons are ideal probes of the initial stages of the 
phase transition as they do not undergo final state interactions and basically freeze out 
with long mean free paths compared to strong interaction length scales. Our strategy will 
therefore be to develop and implement, from first principles, a scheme for calculating the 
number of photons at a given time t after the quench. 

The fact that the photoproduction rates are related to current-current correlators is a 
well-known result. However, in this section, we will derive such a relationship from first 
principles for photoproduction from a bath of pions which is far from equilibrium. 

As argued in the previous section, the lowest order contribution in a arises from the 
non-equilibrium current vertex. Therefore, to this order we neglect both the “seagull” term 
and the Coulomb interaction. We must emphasize that we are performing a perturbative 
expansion in the electromagnetic coupling a , while the pion dynamics is being treated non- 
perturbatively, and in fact exactly to leading order in the 1/IV expansion. Thus to this order 
in a, 


M‘» = -T q i 


(<h T (g) • Jxi—q)) + (Mq) ■ $ T (—q}) 


(4.1) 


The expectation values in the above expression are expressed in terms of correlators of 
operators living on the forward and backward time contours, yielding [fUl] 


(m)) = 


2 ie 2 d 


dt! 


d 3 k 


tfk? 


( 7 T + (k, t)7T + ( — k, t'))x 


q dt" J J (27r) 3 

x (t t + (-k - q, t)k + (k + q, t))($p(-q, t')$T % (q, t")) 

-(tt + (£, t)n ~(-£, t'))(TX + (-k - q, t)n~(k + q, t'))($f J (-q, t')$f l {q, t")) 


(4.2) 


J t=t" 


where we have omitted the isospin indices on the pions and used the fact that (ttiTTi) = (7r 2 7r 2 ). 
The two point functions (7f + (fc, t)n~(—k, t')) = —iG^(t,t') and (7T + (k,t)^ + (~k,t')) — 
—iGf + (t, t') can be read off from (|2.20|) , (|2.24j) and ( [2.25| ). 

Gathering all these ingredients together, we find after some lengthy but straightforward 
algebra 
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(N q <t)) = 


47r 3 g 


dt' / d 3 k[k 2 — (k ■ q) 2 ] x 


(4.3) 


X 


G>(t, t')G> Uql (t, t')g<(f , t) - G<(t, 0 W*) 0(t - 0- 


The theta function ensures that this expressions is causal. The Green’s functions for the 
pions ( |2.24| ), (|2.25| ) are expressed in terms of the non-perturbative mode functions which are 
obtained numerically by evolution of the equations (|2.16|) . Finally the rate can be written 
in terms of the pion mode functions as (see equations ( |2.16l , [2.24l,[2.25|) 


(N q <t)) = 


d 3 k 


-k 2 shr 9 


dt’ 


32vr 3 qJ W k W^ 

{[c/ t («)e;(i')C/| t+ , 1 (i)qE +sl (*')e- ,,( ‘-‘' ) (l + n k )(l + n |E+fl )(l + JV,(0)) 


+2 




|fc+q|^ 


(4.4) 


,(()(/*j +51 ((')e i «<‘-‘ + «| M )JV,(0) 


+ [e;(()f/ l (i')f/ |E+ , 1 («)qj + , 1 («')e-‘«<‘-‘'>n lIl|J+l1 (i + jv,(o)) 

-C/ t (()f/;((')(7 |l+5l («)qE + » 1 (*')e‘ ,( ‘-‘ ,| (l + »u.)(l + nu+flJJV^O)] + c.c.} 

There are several noteworthy aspects of this equation. First, we emphasize that to lowest 
order in perturbation theory in ck, this expression is exact in the specific sense that we have 
not coarse-grained over any time scales in the problem. All the relevant microscopic time- 
scales of the theory are accounted for in this expression. This is a generalized quantum 
kinetic equation, which avoids the assumption of completed successive scattering processes 
implicit in a simple Boltzmann description. There is no energy conservation on short time 
scales (~ a ferrni), but memory of the earlier states remains. This lack of energy conservation 
on short time scales allows one to study “transient phenomena” |39|j4^] . 

Next, the individual terms of this expression have a very physical and simple interpreta¬ 
tion in terms of “rate balancing processes”. As shown in Fig.(4.a) the first term represents 
the difference in rates for two processes: the emission and absorption of a photon and two 
“off-shell” pions in the medium. By “off-shell” pions, we are referring to the pions with 
fully dressed non-equilibrium propagators in the large N approximation. Thus, although the 
processes conserve momentum, the absence of equilibrium means that they are not energy- 
conserving. The second term, depicted in Fig. (4b), represents a “bremsstrahlung” type 
process and the inverse process of “scattering” of a charged pion off a photon. The reason 
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that we identify the direct process with a “bremsstrahlung” type process is because the full 
dressed pion propagator has insertions of the mean-field and the self-consistent fluctuation 
as shown in the figures. These processes can also be identified with a non-perturbative re¬ 
summation (large N) of Landau-damping type contributions, which in equilibrium and in 
the usual perturbative expansion only contribute to processes with four momenta below the 
light cone. Finally the third term, depicted in Fig. (4c) represents 7r + 7r~ annihilation into a 
photon, minus the reverse process, with a photon creating a charged pion pair. 

Another important point is that, as we show below, the rate vanishes automatically if 
the system is in equilibrium: 

In equilibrium, the pion mode functions are positive frequency plane waves i.e. U k (t) = 
exp(— itUkt) = exp(— i\Jk 2 + m 2 t). The time integrals can then be easily performed by intro¬ 
ducing an adiabatic “switch-on” convergence factor for the lower limit of the time integral 
yielding energy conserving delta functions. In particular, we obtain 

• r d ^ k“ c 

mt)) W t 81112 0 + + « )x 

x[(l + n k )( 1 + n^ +?1 )( 1 + N q (Q)) - n k n^N q ( 0)] 


+25(uJk - ^\k+q\ + ?)[(! + n fc) n |fc+5|(l + -Ng(O)) _ n k( 1 + n \k + q])N q {0)} 


+5(u k + u^ -q)[n k n$ + %(1 + N q (0)) - (1 + n k ){\ + n |jf+?| )iV g (0)]} . (4.5) 


The delta functions can never be satished and so the rate simply vanishes in equilibrium, 
due to the kinematics. Furthermore, upon linearizing near the equilibrium distribution for 
the photons IV g (0) = N eqq + SN q , i.e. in the “relaxation time approximation”, one can 
easily see that the expression that multiplies SN q is identihed with the imaginary part of the 
photon self-energy to lowest order in a evaluated on the photon mass-shell, which vanishes 
kinematically. 

This is of course a well known result: the equilibrium photoproduction rate is given by the 
imaginary part of the photon self-energy evaluated on the photon mass shell. In equilibrium 
the lowest order contribution to the on-shell imaginary part of the photon self-energy arises 
from Compton scattering and pair annihilation diagrams (and their reverse processes) in the 
medium; these are both 0(a 2 ) effects. 

Thus we see that out of equilibrium there are “off-shell” processes that give a non-zero 


0(a) contribution to photoproduction. The equation ( |4.4|) and these conclusions are some 
of the more important results of this article. 

Having established the equivalence with the usual result in the case in which the state is 
prepared in equilibrium in the infinite past (in-state) and evolved in time, we now consider 
the kinetics of photoproduction as an initial condition problem. In the situation under 
consideration, i.e. a sudden “quench” below the spinodal region in the linear sigma model, 
the conclusion of the discussion above implies that if one assumes that the system was in 
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(local) thermodynamic equilibrium for all times before the quench t < 0, then the rate 
vanishes for t < 0. Therefore, let us consider the following initial value problem where 
the initial density matrix, prepared at time t = 0 commutes with the respective number 
operators with some initial distribution of photons and let this initial density matrix evolve 
in time with the total Hamiltonian. 

The calculation of the rate goes through in exactly the same manner but now the lower 
limit of the retarded time integral in (|4.4|) is the initial time t = 0 (see appendix A for a 
brief discussion of kinetics as an initial value problem and [|39l,|42|,|43| for more details). 


A. Improved Rate Equation 

In this section we present an improvement to equation ( |4.4| ) obtained in the previous 
section, that represents a resummation of the perturbative series in a. The diagrammatic 
structure of the perturbative series for the rate is recognized to lead to the following exact 
expression for the rate 

(N q <t)) = -£^J q dt ' j - (£ ■ <?) 2 ] x 
K pm (* > (*'> *) - K pm (* > oC (*'> *) 

where Q > ’ < are the full photon Green’s functions, and S >,< are the full photon self-energies 
obtained to all orders in perturbation theory in the electromagnetic coupling. The expres¬ 
sion for the rate given by (|4.4| ) is obtained by replacing Q ’ , S >,< by the zeroth order 

contributions in terms of free held propagators. The exact expression for the rate (|4.6| ) gives 
the rate as a function of time and the initial population N q ( 0), which combined with the 
Schwinger Dyson equation for the photon propagators gives the exact rate. The lowest order 
term, however, neglects the change in the population providing only the lowest order time 
evolution. A resummation must be invoked to incorporate the change of population as a 
function of time. 

The first point to notice is that N q in (|4.4| ) is proportional to the fine structure constant 
a. The presence of this weak coupling is a signal that there is a separation of time scales 
in the problem. To be specific, the time scales we are referring to are: (i) ( N q /N q )~ l which 
is the time scale that governs photoproduction; and (ii) the time scale associated with the 
non-equilibrium processes ~ 1 fm/c which are the time scales for evolution of the pion mode 
functions. The presence of an e 2 in the expression for N q , guarantees that these two time 
scales are widely separated. This means that the photon distribution will show only small 
deviations from equilibrium on short time scales. Obviously these expressions cannot be 
extended to long time scales since they neglect the changes in the initial distribution. 

Now consider implementing the following procedure. Let us integrate (|4.4D from some 
initial time to to a time t 0 + At, (where At < lfm/c, i.e. the microscopic time scale), hence 
obtaining N q (to + At). Assuming that the change in the photon distribution function is small 


0(t — t'). (4.6) 
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during this interval (an assumption warranted for weak coupling), we can then update the 
value of N q that enters into the photon Green’s functions on the RHS of ( [PI) to N q (t 0 + At). 
In terms of the density matrix, this procedure has the interpretation of starting off with a 
density matrix that is diagonal in the number basis. As the system evolves, correlations of 
the form (b q (t)b_ q (t)) and (&£(i)&t 9 (i)) are generated, which give rise to off-diagonal terms in 
the density matrix |39|]. The updating procedure outlined above neglects these correlations 


by collapsing the density matrix to a diagonal one at the end of each infinitesimal time step. 

This procedure is iterated for all times, resulting in the replacement N q (0) —> N q (t'). This 
approximation is similar to that employed in quantum kinetics, known as the generalised 
Kadanoff-Baym approximation J42[]. It can be seen to sum a Dyson-like series for the rate 
43j by writing the formal solution iteratively. 

A similar resummation scheme is implied by the semiclassical Boltzmann equation, in 
which if the occupation numbers are treated in lowest order, the change is linear in time. 
Replacing the occupation numbers by the time dependent ones in the Boltzmann equation 
leads to a resummation and exponentiation of the time series [Ek|. However, as discussed in 


the Boltzmann equation assumes completed collisions that result in a coarse graining in 
time and neglects all of the transient effects and dynamics on short time scales. See appendix 
A and [fh],[4^,[43|] for a more detailed discussion of the approximations involved. 

We note that the pion occupation numbers n k are not updated since the change in the 
occupation numbers for the pions is accounted for by the evolution of the mode functions 
through the Bogoliubov transformation that determines the time evolution of the particle 
number [[19] . 

Since we are using a “box” normalization for the particle states, we now pass to the 
“continuum” normalization (in a volume D) by the replacement 


M) - (2Jr) 


, d«jv(t))/n) 

d 3 q 


(4.7) 


(NM) -> (2tt 


,d((N(t))/(l) _ 3 dR(t) 


= (2tt )‘ 


d 3 q d 3 q 

and obtain the final form of the invariant rate of photoproduction per polarization-. 
- dR(t) a f d 3 k - - rt 


(4.8) 


(27r) 3 1< 7 1—WA = -A— / - k 2 sin 2 6 [ dt! 

1 ] m d 3 q 8tt 2 J W k W lUql J o 

{[C7 t (t)t^(i')£/| SH |(i)^(*')e- <,( ‘-‘' ) (l + 0(1 + «|E + ,i)(l + Ay*')) 


vu*)^(*')q J+ , 1 («)t , l f + , 1 («')^' , -‘') IltIl | t+51 jv 8 ( t ') 


+2 


%(i)^(i')q t+ , 1 (i)l r |t + , n (i')e-‘ 5( ‘-‘' ) (l + n»)n |J+a (l + 
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-C/ fc -(i)C / t(*')C / |S +j1 (()t , |lr + , 1 («> i,|, - , ' ) i! l (l + "iShiW') 




+ n*)(l + »| t - + ,-,)JV»((')] + c.c.} , (4,9) 


where 6 is the angle between k and q. This novel expression for the rate is the basic important 
result of this article. The total photon yield at a time t is given by twice (to account for 
both polarizations) the time integral of (|4.9|). 


V. NUMERICAL STUDY AND RESULTS: 

Having established the proper rate equation that accounts for strong off-equilibrium 
effects we now proceed to a numerical evaluation of the rate and the total number of photons 
produced during the time of the transition. 

The numerical study consists of two stages: in the first stage the evolution of the order 
parameter and the mode functions is solved by integrating the equations ( |2.15| , f2T6| , [2. 17] , [2.18]) 
with m and N — 3 corresponding to 3 pion degrees of freedom. We have chosen to 
represent a “quench” from an initial temperature T t = 1.1T C = 220 MeV to zero temperature. 
This choice has no particular physical significance but serves as an illustration of a quench 
scenario not too far above the critical temperature. We have also performed calculations 
with the initial temperature approching the critical temperature from above. 

This stage of the numerical evaluation produces the evolution of the order parameter and 
provides the mode functions £4(f) for all the values of k considered. Consistent with the 
linear sigma model being an effective theory below 1 GeV, we have kept all k wave vectors 
up to this value. As shown in hgure(l) the order parameter reaches the saturation value 
~ f n on time scales ~ 3 fm/c, this time scale signals the last stages of the phase transition. 
Thus we kept all the mode functions up to this time. 

The second stage of the calculation uses the mode functions as input in the numerical 
evaluation of the rate expression ( |4.9|) . The input for the distribution functions is taken to 
be an equilibrium distribution of photons and pions at the initial temperature which is then 
varied from T* = 1.1T C , down to T c . As we perform the calculation of the rate (|4.9| ), we 
also simultaneously integrate the expression to obtain the total number of photons per unit 
volume at a given time (|4.7| ) (multiplied by |(/|). 

Results: 

The results of the numerical evaluation of the rate are clearly displayed in figures (5- 
8). Figure(5.a) shows the total photon yield (i.e. for both polarizations) per unit volume 
( 27 r) 3 |g| 1/1 (i n units of fm -1 ) at time t =3 fm/c vs \q\ (in units of 200 MeV) . We 

clearly see that the distributions of photons produced during the time of the phase transition 
is out of equilibrium. The long wavelength photons, with energies < 100 MeV can be 
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described by a thermal distribution at a temperature T iw m 1.17 T c . The distribution 
for short wavelength photons is basically not modified from the initial distribution and 
merges with the initial thermal distribution at T sw = 1.1 T c . The distribution function 
smoothly interpolates between a thermal distribution at a temperature which is about 6 % 
above the initial temperature at long wavelengths, and another thermal distribution at the 
initial temperature at short wavelengths, with a smooth crossover at energies ~ 100 MeV. 

To better quantify the production of photons, we show in figure (5.b) the difference 
between the distribution at t — 3 fm/c and the initial distribution, thus subtracting the 
thermal background. The enhancement at low momenta is clearly displayed in this figure 
where it can be seen that the effect is definitely more marked for k < 100 MeV. Figure ( 6 ) 
shows the invariant rate (2Tr) 3 \q\dR(t)/d 3 q (in units of frn -2 ) at time t — 3 fm/c vs |g| (in 
units of 200 MeV). Again, clearly the rate is enhanced at very low momenta, in the same 
range as shown in Figure (5.b) thus explaining the enhancement in the total photon yield 
in this energy range. The reason for this dramatic effect at long-wavelengths is physically 
clear to understand. The pion mode functions that enter in the expression of the rate (|4.9| ) 
grow exponentially because of the spinodal instabilities for k < 100 MeV as shown in figure 
(2.a). Large transferred photon momentum q takes the mode functions outside the band of 
spinodally unstable modes with a much smaller contribution to the rate and total photon 
yield. This phenomenon is quite independent of the initial distribution function for the 
photons and solely a feature of the instabilities associated with the chiral phase transition. 

To see that this result is independent of the initial distribution of photons we show in 
figure (7) (27r) 3 |g| a ^. ^j me 3 f m / c f or the case of initial “vacuum” conditions on the 

photon occupation N q ( 0) = 0. This figure is very similar to the subtracted result displayed 
in figure (5.b), thus explicitly showing that the enhancement is solely a consequence of the 
pion instabilities and independent of the initial photon distribution. Figure ( 8 ) shows the 
invariant rate for this initial condition. 

A further enhancement of the photoproduction rate and total photon yield is obtained 
with quenches closer to the critical temperature, with a dramatic enhancement of about 
15 — 20% at low momentum when the quench is at the critical temperature. This is an 
understandable result because at the critical temperature pions are effectively massless with 
large contribution from their equilibrium distribution functions at long wavelengths. How¬ 
ever, we find that a scenario based on a critical quench is physically quite unlikely. 

At this stage it is very illuminating to compare our results with those obtained for 
photoproduction from a hadronic gas by Kapusta, Lichard and Seibert |I0| and more recently 
by Steele, Yamagishi and Zahed |3T|]. In reference 
photoproduction from qq —> 'yg ; q(q)g —> q(q)g 7 is given by 


30 the result for the invariant rate of 


dR 207T 
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2.9 E 
4na s T 


+ 1 


(5.1) 


which for E < 100 Mev and T ps 200 Mev is at least an order of magnitude smaller 
than the results shown in figures ( 6 , 8 ). I 11 figure 5 of reference |3T| the invariant rate for 
photoproduction at a temperature 150 MeV is shown. Although their initial temperature is 
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smaller than the value used by us, we can see that after normalization of scales and units, 
the rate that we find for momenta q < 100 MeV (shown in figures (6,8)) is several orders 
of magnitude larger than that displayed in that reference. Thus the main conclusion of this 
comparison: the long wavelength instabilities associated with the fast phase transition below 
the critical temperature are responsible for a dramatic enhancement of the photoproduction 
rate and yield for low momenta. 


VI. CONTRIBUTION FROM THE ANOMALOUS DECAY vr° -*• 2 7 : 

One of the most important sources of low energy photons in a typical heavy ion collision 
is the decay of the neutral pion into two photons, with a branching ratio of almost 99%. 
This process typically produces photons with energies >70 MeV. 

Thus it is important that we quantify the contribution of n° decay to the non-equilibrium 
production of photons during the phase transition. It is straightforward to perform the 
calculation leading to the rate of photoproduction, now including the anomalous vertex 
given by ( fi.2| ) . The lowest order contribution is of 0(a 2 ), because each anomalous vertex 
is of 0(a), and depicted in figure (9). Notice that the intermediate state has one pion 
propagator, therefore two mode functions Uk(t). This contribution must be compared to the 
one-loop contribution depicted in figure (3), which is of O(a) and the intermediate state 
has two pion propagators with four mode functions. The important point to notice is that 
the exponentially growing (spinodally unstable) modes will give a much larger contribution 
to the diagram with the pion loop than to the diagram for no decay. Furthermore there is 
a factor of a difference between the two giving an even smaller weight to the neutral pion 
decay diagram in favor of the one-loop pion diagram. 

Thus we conclude that the non-equilibrium process described in the previous section gives 
a far larger contribution to the photoproduction rate during the time of the phase transition 
than the contribution from neutral pion decay. Both because of the fast non-equilibrium 
growth of spinodally unstable modes and because of the powers of a. 

We are then led unequivocally to the conclusion that the non-equilibrium photoproduc¬ 
tion process through spinodal instabilities is far more efficient than 7T° —> 2y for photons 
with energy below 80 MeV. 


VII. DISCUSSION, CONCLUSIONS AND OUTLOOK: 

In this article we have focused on the description of the process of photoproduction during 
the non-equilibrium stages of the chiral phase transition and formation of chiral condensates. 
The premise of this work is that if the chiral phase transition occurs far from equilibrium, 
resulting in the possible formation and relaxation of disoriented chiral condensates, the long 
wavelength pion instabilities will lead to an enhancement of photoproduction at low energies. 

We have developed a novel quantum kinetic approach to photoproduction that accounts 
for the non-equilibrium dynamics in short time scales. This approach incorporates consis- 
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tently the dynamics of long-wavelength pion fluctuations that undergo spinodal instabilities 
during the phase transition. These instabilities lead to an enhanced photoproduction rate 
and photoproduction yield at low energies \q\ <80 MeV. Comparing our results with recent 
results of photoproduction in the hadronic gas lfiO|. fIT| , we are led to conclude that if the 


chiral phase transition occurs far off equilibrium there will be a dramatic enhancement in 
photoproduction in the energy range below 80-100 MeV. This result is independent of the 
initial photon distribution, and its origin resides solely on the strong instabilities of the pion 
fluctuations, and thus on the dynamics of the chiral phase transition. 

Although we have focused on the “quench” scenario, our quantum kinetic approach is 
certainly more general and can be generalized to contemplate the case in which the transi¬ 
tion occurs by bubble nucleation as proposed recently |24|]. Our approach did not consider 
either longitudinal or spherical expansion, since we were interested in understanding if new 
phenomena could emerge. The next step is to extend these methods to include expansion 
|T8| as well as inhomogeneous configurations. In particular, one can now study the fluctua¬ 
tions around semiclassical configurations, such as the Blaizot-Krzywicki DCC and perform 
a calculation along the lines detailed in this article that would include both the classical 
currents and the quantum fluctuations, thus extending the results of [^5| to include the full 
quantum and thermal evolution. 

Another avenue to pursue is a detailed computation of dilepton rates, by extending the 
kinetic approach developed here, also incorporating the effects of hydrodynamic expansion. 
The ultimate goal of such a program is to offer detailed probes of the dynamics of the 
chiral phase transition either by revealing the formation and relaxation of Disoriented Chiral 
Condensates or some other type of non-equilibrium phenomena associated with the chiral 
phase transition. The energy range in which the enhancement occurs is a difficult one for 
the present detectors at AGS and SPS. However both the full acceptance detector (FELIX) 
proposed at CERN-LHC and the PHENIX detector scheduled to operate at RHIC towards 
the end of the millenium may provide the necessary resolution at low energies to provide a 
window to probe these phenomena. 
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APPENDIX A: FORMAL KINETIC EQUATION 


Our treatment of photoproduction far off equilibrium is based on a description of kinetics 
as an initial condition problem. Such an approach, although non-standard in field theory at 
finite temperatures, has become a standard in the studies of fast processes in semiconductor 
physics |42j]. Far off equilibrium situations result in solid state physics when semiconductors 
(or other materials) are studied with fast pulsed lasers, which probe the dynamics of these 
systems on femtosecond scales. It has been recognized in the last few years that a Boltzmann 
approach that coarse-grains over the microscopic time scales misses most of the important 
strongly out of equilibrium effects associated with virtual transitions so that the full quan¬ 
tum kinetic equations must be studied. These equations are typically non-Markovian (with 
memory kernels) and non-linear, to allow for non-linear relaxation, thus going well beyond 
the relaxation time approximation usually invoked in near-equilibrium situations. 

A formal kinetic equation for an initial condition problem is obtained as follows. Consider 
that at some initial time t = t Q the system is described by a given density matrix p(t 0 ) that 
commutes with the number operator N , but evolves in time with a Hamiltonian that can 
depend explicitly on time (as is the case under consideration in this article). The time 
dependent number operator in the Heisenberg picture obeys 


N(t)=i[H(t),N(t)] 


(Al) 


with solution 


N(t) = N(t 0 ) + i f [H(t '), N{t')] dtf 

Jto 


(A2) 


This equation leads to an iterative series for N(t) in terms of higher order commutators 
and N(t 0 ). Iterating this equation once, and using the assumption that the initial density 
matrix commutes with the number operator at the initial time, one obtains the following 
exact expression for the expectation value of the rate in the initial density matrix, 


W)) = ~ ([H(t\[H(t’),N(t')}})dt' 

Jt 0 


(A3) 


The lowest order contribution to the rate results by replacing N(t') by N(t 0 ). The 
expression (|A3|) reveals two important features of the exact expectation value of the rate: 
i) the rate vanishes at t = t 0 as a result of the assumption that the initial density matrix 
commutes with the number operator at the initial time, ii) the number operator enters 
formally in the expression above with the time argument that is being integrated. This 
expression reveals a non-Markovian structure in the rate equation in terms of the memory 
kernels resulting from the nested commutators, and the argument of the operator number is 
always integrated in the retarded time integrals. This observation leads to the generalized 
Kadanoff Bayrn approximation |42| , in which the calculation of the rate is performed to first 
order in perturbation theory, corresponding to replacing N(t r ) = N(t 0 ) in (|A3|) , recognizing 
the proper kernel from this expression (the double nested commutator) and finally replacing 
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the expectation values of the number operator at the initial time by the expectation values at 
the integrated times t'. By expanding the expectation values in terms of number eigenstates 
and inserting the identity in terms of these states one finds that this approximation neglects 
higher order correlations, in particular off-diagonal matrix elements of the time dependent 
number operator in the basis of the number operator at the initial time as discussed in 
section IV. By iterating the resulting expression one finds that this approximation results 
in a resummation of the terms akin to the Dyson series. In fact in equilibrium it can 
be shown [03) that such a resummation is precisely a Dyson-type approximation. Such a 


resummation is also implied in the usual Boltzmann equation, however, in order to obtain the 
Boltzmann equation more drastic approximations must be made. They correspond to setting 
N{t') pc N(t ), taking the contribution from the N(t ) outside the time integral and performing 
the time integral up to t —» oo enforcing energy conservation a la Fermi’s Golden Rule. Such 
a coarse graining approximation is refered to as “the completed collision approximation” 
in the language of quantum kinetics. This approximation completely neglects dynamics on 
short time scales and transient effects. It is known to fail when strongly out of equilibrium 
processes, ocurring on short time scales are important ||42|| . 

Keeping the memory kernels results in a partial resummation of the exact rate equation 
akin to the resummation of particular diagrams in the Dyson series for self energies and 
allows to study “transient” phenomena associated with the initial stages of relaxation of the 
initial state (43j which are completely missed in a Boltzmann description, which leads always 
to exponential relaxation through energy conserving processes. 

Two numerical strategies: 

The typical rate equations of quantum kinetics obtained in this article ( [f.9| ) are non- 
Markovian and therefore non-local in time. There are two ways to deal numerically with 
these : i) carry out the retarded integral directly performing all the momentum integrations 
or ii) recognize that the rate can be written in the generalized form of a sum of factorized 
kernels 


N= dt'J2'Yk(t)K k (t') = r y k (t)H k (t) 


(A4) 


with the supplementary variables H k {t) that obey the differential equations 

= K k {t) ; H k ( 0) = 0 (A5) 

at 

Thus the non-local Markovian kernel is traded for a set of local first order differential equa¬ 
tions. 


APPENDIX B: A PEDAGOGICAL EXERCISE: THE FORCED HARMONIC 

OSCILLATOR 

Since the treatment of kinetics as an initial value problem is not part of the standard lore, 
we present in this appendix a pedagogical exercise with the purpose of providing the reader 
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a simpler setting within which to understand these concepts and that captures the essence 
of the scheme that is used in the field theory calculation. The system that we consider 
is a simple harmonic oscillator coupled to a classical source. This is already a significant 
departure from the actual problem, where the electromagnetic field couples to the charged 
pion current which is a quantum-mechanical object, but will allow us to make contact with 
standard resuls. Thus we proceed with our example which is introduced purely for illustrative 
purposes. The Hamiltonian for our system is 


H 


p 2 uj 2 q 2 

~2 + 2 


+ 3 (t)q 


(Bl) 


where q and p are the canonical coordinate and momenta respectively, while uj is the oscillator 
frequency. The number operator for the oscillator quanta is given by 


N 


1 

u 



uj 2 q 2 


2 


1 

2 


(B2) 


Therefore, 


dN 

dt 


1 

2c jJ 


pp + pp + l xj 2 qq + uj 2 qq 


Using the Heisenberg equations of motion, namely 


q = p, 

P= -u 2 q-j(t) 


we obtain, 


* = -dhp 


UJ 

3(t) 

UJ 


Q- 


(B3) 

(B4) 

(B5) 

(B6) 


^Frorn the equations of motion ([S5j) we see that the Heisenberg operator q satisfies the 
equation 


q + uj 2 q = -j(t). (B7) 

The solutions to the above equation can be written in terms of a homogeneous solution to 
the operator equation and a c-number piece 

q(t) = qo(t) + j G ret (t - t')j(t')dt' (B8) 

where G ret is the retarded Green’s function 
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G ret (t — t') = — sin[o;(i — t')]0(t — t'). (B9) 

UJ 

We now compute operator expectation values in an initial equilibrium state specified at 
t = t 0 which we chose to be the ground state of the unperturbed harmonic oscillator in 
which (p) = (q) = O.The current is switched on at time t = t 0 . Considering this initial state, 
the expectation value of the number operator in the time evolved state in the presence of 
the current is given by 

(ft) = — ^-^-(q) = — - [ dt 1 cos[u;(f — t')]j(t') (BIO) 

UJ UJ Jto 

leading to the final result for the number of quanta produced up to time t 

{N(t)) = —[ dti [ dt 2 j(ti)cos[uj(t 1 -t 2 )]j(t 2 ) (Bll) 

UJ J to J tQ 

We can evaluate the same expectation value in the closed time path formalism outlined in 
the previous section i.e. 


(ft) = f V[q + }V[q-}q + (t)eI dt ^ L+ - L ^ 

l o dt J 

where the lagrangian is 

L = ^q 2 - i u 2 q 2 -j(t)q. 

Therefore, 

(ft) = _MA( g +(t) e -*/^(t')(9 + (C-?-(C)) 0 . 

uj dt 


(B12) 


(B13) 


(B14) 


Expanding out the exponential and imposing the tadpole condition i.e. (q + ) 0 = 0, we see 
that the only non-zero contribution is from the first-order term, which yields: 


(ft) = i — 4 / dt , j(t')((q + (t)q + (t ')) 0 - ( q + (t)q (t') 0 ). 

uj dt J to 

The two point functions can be read off from ( 2.2C ), ( 2.22|) and ( 2.24 ), giving 



f dt'j(t')(e- ibj{ t-^ 

2 u 2 dt Jt 0 

/ dt'jUjcoslujU — t 1 )]. 

UJ J to 


(B15) 


(B16) 


Thus first order perturbation theory gives us the exact result. This is of course due to the 
fact that the current is a c-number object. 














At this point we can make contact with the usual in-out approach based on S-matrix 
theory to calculate rates. Consider the expression (|B11| ) in the limit to —> — oo, the total 
number of quanta created at t —> oo is given by 


(iV(oo)) = 


\Ml 

2u 


(B17) 


with j(cj) being the Fourier transform of the current evaluated at the oscillator frequency 
u. This is the standard result for the case of a non-interacting theory in the presence of a 
c-number current and makes the connection with the treatment of Huang and Wang (2^j for 
the sigma model. 
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Figure 3 


FIG. 3. Diagrams contributing to 0(a) to the non-equilibrium photoproduction rate. 
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Figure 4(a) 



Figure 4(b) 



Figure 4(c) 

FIG. 4. Processes that contribute to the photoproduction rate. The pion propagators are the 
full large N-resummed propagators. 
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FIG. 7. Total photon yield per unit volume (27r) 3 |g|in units of fm 1 vs. |g| in units 
of 200 Mev, at time t =3fm/c, for “vacuum” initial photon occupation, N q (0) — 0. 
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FIG. 8. The invariant photoproduction rate (27 t) 3 |< 7 | d in units of fm 2 vs. 
Mev at time 3 fnr/c, for “vacuum” initial photon occupation, N q ( 0) = 0. 
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Figure 9 


FIG. 9. Lowest order contribution to photoproduction from the anomalous decay of the neutral 
pion via Ca in (3.2). 
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